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ABSTRACT 

The field containing the candidate High Mass X-ray Binary IGR JO 1363+66 10 was observed by XMM-Newton 
on 2009 July 31 for 28 ks. A Be star was previously suggested as the possible counterpart of the INTEGRAL 
source, and although Chandra, during a 2007 observation, did not detect an X-ray source at the position of 
the Be star, we find a variable source (XMMU J013549.5+661243) with an average X-ray flux of 2 x 10"° 
ergs cm -2 s" 1 (0.2-12 keV, unabsorbed) at this position with XMM-Newton. The spectrum of this source is 
consistent with a hard power-law with a photon index of Y = 1.4 ±0.3 and a column density of A^h = (LS^'J) x 
10 22 cm" 2 (90% confidence errors). These results, along with our optical investigation of other X-ray sources 
in the field, makes the association with the Be star very likely, and the 2 kpc distance estimate for the Be star 
indicates an X-ray luminosity of 9.1 x 10 31 ergs s" 1 . This is lower than typical for a Be X-ray binary, and 
the upper limit on the luminosity was even lower (<1.4 x 10 31 ergs s" 1 assuming the same spectral model) 
during the Chandra observation. We discuss possible implications of the very low quiescent luminosity for the 
physical properties of IGR JO 1363+66 10. 

Subject headings: stars: neutron — X-rays: stars — stars: emission-line, Be — black hole physics — stars: 
individual (IGR J01363+6610) 



1. INTRODUCTION 

The hard X-ray imaging of the Galactic plane by the 
International Ga mma-Ray Astrophys ics Laboratory (INTE- 
GRAL) satellite dWinkler et al.ll2003l) has uncovered a large 
number of new or previously poorly st udied "IGR" sources 
dBodaehee et all 120071 iBlrJetan 12010). While INTEGRAL 
excels at detecting sources in the 20-50 keV band, it only 
localizes the sources to l'-5', requiring follow-up obser- 
vations with other X-ray satellites to obtain secure optical 
or IR counterparts, allowing for a det e rmination of the na- 
ture of the sources dWalter et all 120061: iTomsick et al.1 120081 
iRodriguez. Tomsick & Chatyll2009l) ~ 

As more and more of these sources have been identi- 
fied, possibly the biggest surprise is the large number of 
High Mass X-ray Binaries (HMXBs) as well as the proper- 
ties of the se systems. Many o f the dozens of INTEGRAL 
HMXBs dBodaghee et al.1 [20071) have large lev els of intrin- 
sic absorption with A^h ~ 10 23 -10 24 cm" 2 (e.g.. IWalter et al.l 
2006), and these are commonly called obscured HMXBs. In 
many cases, it appears that this is due to th e compact object 
being embedded in a strong stellar wind (iFilliatre & Chatvl 
12004b iMoon et al.ll2007t IChaty et al.ll2008l) . Some members 
of the group of INTEGRAL HMXBs exhibit other extreme 
properties, including the high amplitude X- ray flaring of the 
Supergiant Fast X -ray Transients (SFXTs, I in' t Zand! 120051: 
Smith et al. 2006) or long-period pulsations from very slowly 
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rotating neutron stars dPatel et al.ll2007l) . 

INTEGRAL observations have also led to the addi- 
tion of more HMXBs in the B e X-ray binary class 
dRappaport & van den Heuvellll982l) . The optical flux from 
these systems is dominated by an early-type star with emis- 
sion lines from a circumstellar disk. In most cases, tran- 
sient X-ray emission demonstrates the binary nature of the 
system as eccentric orbits lead to periodic X-ray outbursts 
when the compact object approaches periastron. Of the 
64 known Be X-ray binary systems, X-ray pulsations in- 
dicate the presence of a neutron star in 42 cases, and the 
compact object type is unknow n for the remaining systems 
dBelczynski & Ziolkowskil2009t) . One of the interesting prop- 
erties of this class is the observed correlation between the or- 
bital period and the spin period of the neutron star dCorbetl 
1986). 

IGR JO 1363+66 10 has been tentatively classified as a Be 
X-ray binary. The source was discovered during INTEGRAL 
observations on 2004 April 19 but was not det ected ~2 weeks 
later, indicating that the source is transient (Greben ev et al.1 
2004). The peak flux observed from the source was 17 mcrab 
(~2.6 x 10" 10 ergs cm" 2 s" 1 ) in the 17-45 keV band and 9 
mcrab (~9.1 x IP" 11 ergs cm" 2 s" 1 ) in the 8-15 keV band 
dGrebenev et al.ll2004l) . The 3 -a upper limit in the 17—45 keV 
band 2 weeks later was <11 mcrab (Grebe nev et al.l 120041) . 
and the source has not been detected in other INTEGRAL ob- 
servations even though 2.3 Ms of INTEGRA L time have bee n 
accumulated at the position of this source dBird et al.ll2010h . 
Reports of non-det ection include both d etailed studies of the 
Cassiopeia region (Iden Hartog et al.l[2006l) and catalogs indi- 
cating that the source was only detected during the discov- 
ery o utburst in 2004 dKrivonos et al.1 120071: iBird et al.l 120071 
l2010h . 

Although the INTEGRAL position uncertainty of 3'.7 leaves 
a large error region, a Be star was found within the error 
circle using narrow-band Ha imaging and follow-up opti- 
cal spectroscopy, and it has been suggested to be the likely 
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counterpart dReig et a l. 2004, 2005). However, a sensitive X- 
ray observation taken with the C handra X-ray Observ atory in 
2007 failed to detect the Be star dTomsick et al.ll2008l). A t the 
2 kpc distance estimated for the Be star (Reig et al. 20Q5|), the 
non-detection implies an uppe r limit on the X-ra y luminos- 
ity of <2 x 10 31 ergs cm" 2 s" 1 dTom sick et al. 2008), which is 
lower than quiescent luminosities for other Be X-ray binaries 
(Cam pana et alj|2002l) . This luminosity approaches the level 
that has been seen during quiescent periods from transient 
Low Mass X-ray Binaries (LM XBs) due to thermal emission 
from t he neutron star surfac e jBrown. Bildsten & Rutledge 
11998b ICampana etail 119981) . Although neutron stars in 
HMXBs might not be heated to the high levels seen for 
LMXBs, it is interesting that these observations probe this lu- 
minosity regime. 

Currently, there are significant uncertainties about the na- 
ture of IGR JO 1363+66 10. While we know that it is an X-ray 
transient, and the hard X-ray emission makes it likely that it is 
a binary (although its orbital period is unknown), the Chandra 
non-detection makes it unclear whether it is really a Be X-ray 
binary. Finding a Be star in the relatively large INTEGRAL 
error circle is not convincing because Be stars are more com- 
monly f ound as single stars rath er than being part of a binary 
system dPorter & Riviniusll2003l) . Furthermore, we cannot be 
certain that the compact object in the system is a neutron star 
since pulsations have not been detected. 

In this paper, we report on a second sensitive X-ray ob- 
servation of the IGR JO 1363+66 10 field with XMM-Newton 
along with optical spectroscopy of the Be star as well as other 
X-ray sources with optical counterparts in the field. With 
XMM-Newton, we confirm that the Be star is an X-ray source. 
We also re-analyze the Chandra observation and discuss the 
2007 results in the context of the new information from XMM- 
Newton. 

2. OBSERVATIONS 

2.1. XMM-Newton 

We observed the IGR JO 1363+66 10 field with XMM- 
Newton on 2009 July 31 from 14.4 h to 22.3 h UT. The obser- 
vation (ObsID 0603850101) occurred during XMM-Newton 
revolut ion 1766. The EPIC pn, MOS1 and MOS2 instru- 
ments dStruderet all 120011: iTurneret af] 1200 ll) were all op- 
erated in Full Frame mode with a medium filter. We used 
the XMM-Newton Science Analysis Software (SAS-10.0.0) 
package to process the raw data files. We used the SAS tool 
cif build to obtain the necessary "current calibration files" 
for the observation. Then, we reprocessed the pn, MOS1, and 
MOS2 data using epproc and emproc, yielding photon 
event lists for the three instruments. Although proton flares 
sometimes cause portions of XMM-Newton observations to 
have very high backgrounds, we did not find evidence for pro- 
ton flares, and we were able to use the full exposure time. 

2.2. Chandra X-ray Observatory 

We also used an observation of the IGR JO 1363+66 10 field 
made by the Chandra X-ray Observatory. The observation 
(ObsID 7533), which had an exposure time of 4,976 s, was 
made on 2007 June 8 . We used the Adv anced CCD Imaging 
Spectrometer (ACIS; Garmire et al. 2003), and the aimpoint 
was placed on the ACIS-S3 chip. The 90 % confidence IN- 
TEGRAL error circle for IGR JO 1363+66 10 dBird et al.ll2010l) 
is contained on the ACIS-S3 and ACIS-S4 chips. Results 
of the observation were previously reported in Tomsi ck~et al.l 



(2008), and the procedures used for data reduction are de- 
scribed in that work. 

2.3. MDM Optical Observations 

Optical spectroscopic observations were made with the 
Boiler & Chivens CCD spectrograph on the 2.4 m Hiltner 
Telescope of the MDM Observatory on 2009 August 23 UT 
Conditions were clear, but with high humidity. The 150 g/mm 
grating, blazed at 4700 A, was used with a 1" wide slit, which 
gave a wavelength coverage of 3700-7360 A at 7.6 A reso- 
lution. Spectra of six stars coinciding with X-ray sources in 
the field of IGR JO 1363+66 10 were obtained, each with an 
exposure time of 120 s. Of these, three are inside the refined 
INTEGRAL error circle and are presented here. Spectral re- 
duction was per formed using standar d IRAF procedures. Flux 
calibration used Oke & Gunn ( 1983) standard stars, although 
the spectrophotometry is not expected to be accurate because 
of the narrow slit. 

3. RESULTS 

3.1. X-ray Sources in the INTEGRAL Error Circle 

For the XMM-Newton analysis, we used the SAS tool 
edetect_chain to search for sources in the ~30' diame- 
ter field-of-view (FOV) of the MOS1 and MOS2 instruments, 
and included photons in the 0.1-10 keV energy band. All 
seven of the MOS2 CCDs were operational, but, for MOS1, 
CCD#6 was not active because it was switched off after an 
anomaly that occurred in 2005. While this means that some 
of the MOS2 FOV is not covered by MOS 1 , the central CCD, 
which includes the entire INTEGRAL error circle, is covered 
by both MOS units. 

We found a total of 21 sources, including 14 detected by 
both MOS units and seven detected by MOS2 in the part 
of the FOV that MOS1 did not cover. We determined the 
number of counts for each source using 25" radius apertures 
and subtracted the background using circular apertures of the 
same size in 16 source-free regions of the detector. The 
mean number of background counts per aperture is 47.1 ± 1 .8 
for MOS1 and 45.9 ± 1.7 for MOS2 (in 28,180 s of expo- 
sure time). The sources range in brightness from 12 ±8 to 
156 ± 13 counts (both of these sources are in the region cov- 
ered only by MOS2). Five of the sources are in the 3'.7 IN- 
TEGRA L 90% confidenc e error circle for IGR JO 1363+66 10 
given in Bird et al.l d2010t) . These five sources include the sec- 
ond brightest source in the field (132 ±9 counts) as well as 
sources with 74 ± 8, 3 1 ± 6, 25 ± 6, and 19 ± 6 counts, where 
these numbers of counts are averages of the two MOS detec- 
tors. 

The XMM-Newton source names and positions are given in 
TableQ] In addition to MOS1 and MOS2, we also determined 
the position measured by the pn instrument for each source, 
and we report the average position measured by the three in- 
struments (except for XMMU J013632.4+660924, which fell 
between two pn CCD chips). In averaging the positions, we 
weighted each measurement by its statistical uncertainty. In 
Table Q] we also report the overall 90% confidence position 
uncertainties, and, in each case, the error is dominated by the 
systematic pointing uncertainty of 3".4 6 

6 In the document entitled "EPIC status of calibration and data analysis" 
(XMM-SOC-CAL-TN-0018), Guainazzi et al. report an rms value for the 
systematic pointing uncertainty of 2".0, and we have multiplied this by 1.7 to 
obtain the 90% confidence value. 
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Table 1 

XMM-Newton SOURCES IN THE INTEGRAL ERROR CIRCLE 



Name R.A. (J2000)" Decl. (J2000) Position MOS 

Uncertainty'' counts^ 



XMMU J013549.5+661243 


01 


'35' r 


49 s 


53 


+66° 12 


43" 


1 


±3" 


4 


132±9 


XMMU J013606.5+661304 


01 


, 36 ,r 


06 s 


54 


+66° 13 


04" 


5 


±3" 


9 


25 ±6 


XMMU J013620.8+660851 


01 


'36 w 


20 s 


80 


+66° 08 


51" 





±3" 


7 


31 ±6 


XMMU J013632.4+660924"' 


01 


'36" 


32 s 


48 


+66° 09 


24" 





±4" 


1 


19±6 


XMMU J013644.2+661302 


01 


'36" 


44 s 


26 


+66° 13 


02" 


2 


±3" 


5 


74±8 



"The position is the weighted average of the positions measured by the three XMM-Newton instruments. 

*The 90% confidence uncertainty in the XMM-Newton position. This includes a systematic contribution of 3". 4 due to the absolute pointing uncertainty and a 
statistical contribution. We have added the two contributions in quadrature. 
c The average of the 0.1-10 keV count rates measured by MOS1 and MOS2. 

rf This source is detected by MOS1 and MOS2, but not by the pn instrument. It falls on a gap between pn CCD chips. 



Using the 2007 observation of the field by Chandra, we 
pre viously reported three sources in the INTEGRAL error cir- 
cle dTomsick et al.l 120081). However, in th e most recent IN- 
TEGRAL source catalog (Bird et al. 2010), the best estimate 
of the IGR J01363+661 position has shifted by 1'.8 rela- 
tive to the value used in iTomsick et al.l d2008l) . and a fourth 
Chandra source, CXOU JO 13644. 5+66 1301, is now also con- 
sistent with the INTEGRAL position for IGR JO 1363+66 10, 
and its Chandra position is R.A. (J2000) = 01 h 36 m 44 s .54, 
Decl. (J2000) = +66°13'01".6 with a 90% confidence uncer- 
tainty of 0".64. Of the four Chandra sources, two are coin- 
cident with XMM-Newton sources. Thus, merging the 2007 
Chandra source list with the list of 2009 XMM-Newton leaves 
a list of seven X-ray sources in the INTEGRAL error circle, 
and these sources are listed in Table [2] 

3.2. Optical and IR Counterparts and Optical Spectroscopy 

For the seven X-ray sources listed in order of R.A. in Ta- 
ble|2j we searched for optical and infrared (IR) counterparts in 
the United States Naval Observatory (USNO-B1.0) and 2 Mi- 
cron All-Sky Survey (2MASS) catalogs. Four of the sources 
(designated as #1, #2, #5, and #7) have optical and IR coun- 
terparts, and the names of these counterparts and their magni- 
tudes are given in Table [2] For the other three X-ray sources 
(#3, #4, and #6), the nearest USNO and 2MASS sources are 
>7" away, which indicates that these X-ray sources are not 
associated with any sources in the USNO or 2MASS catalogs. 
FigureQ]shows a red (close to /?-band) optical image from the 
Digitized Sky Survey with the revised INTEGRAL error circle 
and the seven X-ray sources labeled. 

X-ray source #1 is the Be star that was previously iden- 
tified using optic al imaging and spectroscopy taken in 2004 
dReig et alll2004 120051) . The optical spectrum of this star 
(USNO-B1.0 1562-0030282) is shown in Figure [2] and has 
the blue continuum and the Ha and H/3 emission lines in- 
dicative of a Be star. The equivalent width (EW) of Ha is 
-54 ± 3 A , whic h is consistent with the value measured by 
IReig et all d2004l) . suggesting that the Be star's circumstellar 
disk is stable. Interstellar absorption features in the spectrum 
can be used to estimate the extinction E(B-V) following the 
correlations in iHerbia d 19751) . Most commonly, the 4430 A 
diffuse interstellar band is used for this purpose, but it falls 
in a poorly exposed region of our spectrum. Instead we use 
the 5780 A feature. With an EW of 1 .0 A, it corresponds to 
E(B-V) in the range 1 .5-2 .0, which is consistent with the es- 
timate of Reig et al. (2005) based on the spectral classification 
and photometry, E(B - V) = 1.6. 



• N 
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FIG . 1 . — Red optical image from the Digit ized Sky Surv ey. The large circle 
is the 3'.7 INTEGRAL error circle from lBird et alj 120101) . The seven X-ray 
sources detected by XMM-Newton and Chandra are labeled. The arrows point 
to the optical counterparts for four of the X-ray sources, and their optical and 
IR magnitudes are given in Tableff] The small circles mark the locations of 
the X-ray sources without known optical counterparts. 



X-ray source #2 is coincident with USNO-B1.0 1562- 
0030364, which is also known as the bright (V = 11.5) op- 
tical source TYC 4043-860-1. There is some confusion about 
whether this is an emi ssion line star based on a cata log of stars 
with Ha in emission (Gonzalez & Gonzalez 1956). However, 
Ha images and optical spectra taken, respectively, 2 months 
and 5 months after the X-ray outburst from IGR JO 1363+66 10 
do not show evidence for Ha emissio n and the spectrum 
shows Ha in absorption (IReig et alj2004l) . The XMM-Newton 
observation provides the first evidence that this is also an X- 
ray source, but the lack of an Ha emission line in 2004 makes 
it very unlikely to be the correct IGR JO 1363+66 10 counter- 
part. 

We also obtained optical spectra for the only other two 
X-ray sources in the INTEGRAL error circle with opti- 
cal/IR counterparts. The spectra for sources #5 and #7 
are shown in Figure [3] and we identify them as G-type 
and early M-type stars, respectively. Like Be X-ray bina- 
ries, LMXBs usual ly show Ha in emissio n during outbursts 
dFender et al.l2009l) and during quiescence dOrosz et all2002t 
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Table 2 

Optical/Infrared Identifications 



Catalog/Source" 


Separation 6 Magnitudes 






(1) XMMU J013549.5+661243 


USNO-B1.0 1562-0030282 


2".0±3".4 5=14.9 ±0.3 R = 12.4 ±0.3 


/ 


= 1 1 .1 ± 0.3 


2MASS J01354986+6612433 


2".0±3".4 J= 10.04 ±0.02 H = 9.57 ±0.03 


K s 


= 9.12 ±0.02 


(2) XMMU J013606.5+661304 


USNO-B1.0 1562-0030364 


2".5±3".9 B=12.1±0.3 /?= 10.7 ±0.3 


I 


= 10.1 ±0.3 


2MASS JO 1360684+66 13021 


3".0 ± 3".9 J= 10.42 ± 0.02 H = 10. 10 ± 0.03 


K s 


= 10.00 ±0.02 


(3) CXOU J013609.9+661157 


The closest USNO-B1.0 source is 9".l ±0".6 away. 






The closest 2MASS source is 


?".0 ± 0".6 away. 






(4) XMMU JO 13620. 8+660851 


The closest USNO-B 1 .0 source is 7".2 ± 3".7 away. 






The closest 2MASS source is 


B".7±3".7 away. 






(5) CXOU J013621.2+660928 


USNO-B1.0 1561-0031113 


0".52±0".64 B = 17.2 ±0.3 R = 13.6 ±0.3 


I 


= 12.0 ±0.3 


2MASS J013621 15+6609286 


0".38±0".64 J= 10.42 ±0.03 H = 9.54 ±0.03 


K s 


= 9.26 ±0.02 


(6) XMMU J013632.4+660924/CXOU J013632.8+660924 


The closest USNO-B 1 .0 source is 9". 1 ± 0".6 away. 






The closest 2MASS source is 7".5 ± 0".6 away. 






(7) XMMU JO 1 3644.2+66 1302/CXOU J013644.5+661301 


USNO-B 1.0 1562-0030589 


0".73±0".64 B= 17.4 ±0.3 R = 14.4 ±0.3 


I 


= 13.1 ±0.3 


2MASS J01364458+6613014 


0".33±0".64 J= 12.77 ±0.02 H= 12.01 ± 0.03 


K s 


= 11.83±0.03 



"The catalogs are the 2 Micron All-Sky Survey (2MASS) and the United States Naval Observatory (USNO-B 1.0). 

''The uncertainty on the separation is the addition (in quadrature) of the X-ray source position error and the 2MASS or USNO-B 1.0 error. 
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FIG. 2.— Optical spectrum of the Be star XMMU J013549.5+661243 taken 
on 2009 August 23, which is a few weeks after the XMM-Newton observation. 
The optical spectrum was taken with the MDM 2.4 meter telescope and shows 
strong Ha and H/3 lines, indicating the presence of a circumstellar disk. The 
lower spectrum is the upper spectrum divided by 6, allowing for all of the 
Ha line to be visible. This source was previously suggested as the most 
likely counterpart to IGR J01363+6610, and the detection of the Be star with 
XMM-Newton provides further confirmation of this association. 



ICharles & Codl2006l) . Also, symbiotic systems oft en show 
strong Ha in emission (Chakrabarty & Roche1 [T997l) . How- 
ever, there are counter-examples for both LMXBs and symbi- 
otics. In any case, as sources #5 and #7 do not have emission 
lines, there is no reason to consider that they might be the 
IGR JO 1363+66 10 counterpart. 

Thus, these observations show that the Be star that 
was previously considered to be the likely coun- 
terpart of IGR JO 1363+66 10 is coincident with 
XMMU J013549.5+661243, which was the brightest X- 
ray source in the INTEGRAL error circle during the 2009 
XMM-Newton observation. The fact that we have now 
conclusively shown that the Be star is a transient X-ray 
source and that none of the other fainter X-ray sources 
in the INTEGRAL error circle have optical properties that 
would be expected of a hard X-ray source strengthens the 
association between the Be star and IGR JO 1363+66 10, and 




4000 4500 5000 5500 6000 6500 

Wavelength (A) 

FIG. 3. — Optical spectra of two of the X-ray sources in the INTEGRAL 
eiTor circle for IGR JO 1363+66 10. The stars are listed as source #5 (top 
panel) and source #7 (bottom panel) in Table [2] We identify source #5 as a 
highly reddened G-type star and source #7 as an early M-type star. 



we focus exclusively on the X-ray properties of the Be star 
XMMU J013549.5+661243 in the remainder of this paper. 

3.3. XMMU J013549. 5+661243 

3.3.1. Energy Spectrum 

We used the SAS tool xmmselect to produce MOS1, 
MOS2, and pn energy spectra for XMMU J013549.5+661243 
and considered the recommendations given in a recent EPIC 
(MOS and pn) calibration document 7 for event filtering and 
energy ranges. For MOS, we used event filtering with the 

7 The document entitled "EPIC status of calibration and data analysis" by 
Guainazzi et al. (XMM-SOC-CAL-TN-0018) is based on results obtained 



The Be X-ray binary IGR JO 1363+66 10 in quiescence 
(a) XMM/MOS1, 2009, 0.1-10 keV (b) Chandra/ACIS, 2007, 0.3-10 keV 




o 



60 arcsec 



FIG. 4. — (a) X-ray image from the 2009 July 31 XMM-Newton observation showing the detection of XMMU J013549.5+661243, which is the very likely 
counterpart to IGR J01363+6610. The image uses 0.1-10 keV photons from the MOS1 instrument and was rebinned so that each pixel is 4" wide. The inner 
circle is the source region used for making the spectrum and the light curve. The annulus (outer two circles) is the region used to estimate the background, (b) 
X-ray image from a 2007 Chandra observation. The image uses 0.3-10 keV photons from the ACIS instrument and was rebinned so that the pixels are 1". The 
circle shows the 90% confidence XMM-Newton error circle, highlighting the fact that IGR JO 1363+66 10 was not detected. In both images, North is up and East 
is to the left. 



expression "#XMMEA_EM && PATTERN <= 12" and in- 
cluded events within an aperture with a 25" radius. A back- 
ground spectrum was extracted from an annulus centered on 
the source, and the source and background regions are shown 
in Figure [4^. For the pn instrument, we used event filter- 
ing with the expression "FLAG=0 && PATTERN <=4" and 
extracted spectra from source and background regions. For 
all three instruments, we used the SAS tools rmf gen and 
arf gen to make response matrices. For the MOS detectors, 
we used the 0. 1-10 keV bandpass, and for the pn detector, we 
used the 0.2-12 keV bandpass. We rebinned each spectrum, 
and we fitted the spectra using x 2 statistics. 

We used the XSPEC version 12 software for spectral fit- 
ting, and we tried several models. The results of using an ab- 
sorbed power-law are shown in Tableland Figure[5] For ab- 
sorpti on, we used the photoelectric absorption cross sections 
from iBalucinska-Ch urch & McCammon ( 1992) and elemen- 
tal abundances from IWilms. Allen & McCrayl (120001) . which 
correspond to the estimated abundances for the interstellar 
medium. This model gives a good fit (x 2 /V = 23.1/26), 
and requires a relatively hard power-law photon index of 
r = 1.4 ±0.3. The power-law fit gives a 0.2-12 keV un- 
absolved flux of (I.9+Q2) x 10~ 13 ergs cm" 2 s" 1 and a col- 
umn density of Nn = (1.5_Q5) x 10 22 cm" 2 , which is some- 
what higher than the value throug h the Galaxy along th e line 
of sight, N H = 5.2 x 10 21 cm" 2 dKalberla et all [2005). Al- 
though this could indicate a small amount of absorption lo- 
cal to the source, it is much lower than the values near 10 23 
cm" 2 seen for the sources usually considered to be obscured 
HMXBs. Another argument against local absorption of the X- 
ray source is that the value of E(B-V) = 1 .5-2. correspond s 
to an jVh in the range (1.0-1.4) x 10 22 cm" 2 dRvteri [19961) . 

using SAS-10.0.0. It was released on 2010 July 16 and can be found at 
http://xmm2.esac.esa.int/external/xmm_sw_cal/calib/index.shtml 




Energy (keV) 

FIG. 5.— XMM-Newton energy spectrum for XMMU J013549.5+661243 
fitted with an absorbed power-law model. The model parameters are given in 
Table[3] indicating a relatively hard spectrum (T = 1 .4 ± 0.3) but not a highly 
absorbed spectrum. The pn data are shown in black, and the MOS data are 
shown in grey with diamonds plotted on MOS 1 points and squares plotted on 
MOS2 points. 



which is consistent with the column density determined from 
the X-ray spectrum. 

We also fitted the energy spectrum with blackbody and ther- 
mal Bremsstrahlung models, and the parameters are given in 
Table [3] A fit with a 1 .3 keV blackbody is similar in qual- 
ity to the power-law (\ 2 / v = 2\.5j 26), but the spectrum does 
not allow us to distinguish between thermal and non-thermal 
models. However, at 2 kpc, the blackbody model implies an 
emitting region with a radius of only 1.15 x 10 3 cm, which is 
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too small to be physical. The Bremsstrahlung model requires 
a relatively high temperature (> 1 1 keV), which makes it very 
similar to a power-law over the XMM-Newton bandpass. Al- 
though we are not able to distinguish between models, it is 
important to note that the power-law fit above shows that the 
spectrum is hard, as expected for an IGR source and espe- 
cially for an HMXB. 

3.3.2. Long- and Short-Term X-ray Variability 

Comparing the two panels of Figure|4]shows the difference 
between fluxes in 2009 (the XMM-Newton observation) and 
2007 (the Chandra observation). We did not detect the Be 
star d uring a 5 ks Chandra/ AOS observation ( Tomsi ck et al.l 
120081) . and 1 count is detected by ACIS within the XMM- 
Newton error circle. Thus, using Poisson statistics, this cor- 
responds toa 90% confidence upper limit of <3.9 counts 
(Gehrels 1986). For the power-law spectral model mea- 
sured by XMM-Newton (T = 1.4 and Nn = 1.5 X 10 22 cm" 2 ), 
this corresponds to an absorbed 0.3-10 keV flux of <1.6 x 
10~ 14 ergs cm" 2 s" 1 and an unabsorbed flux of <2.9 x 10~ 14 
ergs cm" 2 s" 1 , which is >6.6 times lower than the flux seen 
during the XMM-Newton observation. During the original 
outburst detected by INTEGRAL in 2004, an 8-15 keV flux of 
9 mcrab (9.1 x 10" 11 ergs cm" 2 s" 1 ) was measured using the 
JEM-X instrument. In the 8-15 keV band, the flux measured 
with XMM-Newton in 2009 was 7.9 x 10" 14 ergs cm" 2 s" 1 , 
which is a factor of 1,150 lower than the outburst level. Thus, 
the quiescent upper limit obtained with Chandra represents a 
flux >7,600 times lower than the level seen during outburst. 

We also find evidence for shorter term variability during 
the XMM-Newton observation. Figure [6] shows the source and 
background rates combined for all three XMM-Newton detec- 
tors in the 0.2-12 keV band with 250 s time bins. The light 
curve shows that the 0.2-12 keV rates can change from values 
as low as -0.003 ±0.007 to 0.061 ±0.018. We used a \ 2 test 
to determine the significance of the variability. The best fit 
that can be obtained with a constant rate has a x 2 = 142.5 for 
106 degrees of freedom, which indicates that the variability is 
significant at the 99% confidence level. 

4. DISCUSSION 

4.1. The Nature of IGR J01363+6610 

Our results show that XMMU J013549.5+661243 is a 
hard X-ray source coincident with a Be star at a distance 
of ~2 kpc. At this distance, the inferred X-ray luminos- 
ity measured by XMM-Newton is 9.1 x 10 31 ergs s" 1 . While 
this would be a relatively high X-ray luminosity for an iso- 
lated Bl star, it is not too far above t he values reported by 
ICohen. Cassinelli &^ Macfarlane ( 1997) for such stars. How- 
ever, in nearly all cases of isolated B-type stars, the X- 
ray emission is thermal wit h temperatures that are typically 
not above a few xlO 6 K (ICohen. Cassinelli & Macfarland 
Il997t ICohenll2000t) . and their spectra fall very steeply above 
^0.25 keV. The very hard X-ray spectrum that we measure 
for XMMU J013549. 5+661243 is inconsistent with a thermal 
spectrum with such a low temperature, indicating that the Be 
star must have a compact binary companion that is emitting 
most of the X-ray emission. 

While the X-ray luminosity and spectrum are not consis- 
tent with an isolated star, the XMM-Newton data alone are 
not sufficient to determine the nature of the compact binary 
companion. There are known Be systems with both white 
dwarf and neutron star companions that can produce hard 




10000 15000 20000 25000 
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FIG. 6. — Source (a) and background (b) XMM-Newton light curves of 
XMMU J013549.5+661243. The rates are a combination of counts from all 
three instruments (pn, MOS1, and MOS2), and the background rates have 
been normalized to the rates predicted in the source region. The energy range 
in both panels is 0.2-12 keV, and the time resolution is 250 s. The source 
shows variability on this time scale. 



spectra. For example, there are several 7 Cas-like systems that 
may harbor white dwarfs and have powe r-law spectral com - 
ponents with photon indices of 1.4-1.7 dMotch et al.1 12007). 
However, based on our study of the X-ray sources in the 
IGR JO 1363+66 10 field, we argue that it is very likely that 
XMMU J013549.5+661243 is the quiescent counterpart to 
IGR JO 1363+66 10. In this case, the outburst detected by IN- 
TEGRAL indicates that the source can produce an X-ray lu- 
minosity as high as ~40 35 ergs s" 1 along with variations in 
X-ray luminosity by a factor of >7,600 (when the Chandra 
observation is considered as well). In contrast, the 7 Cas-like 
systems have luminosities typically near 10 33 ergs s" 1 , and 
while they show variability in flux by factors o f a few, they 
do not have large outbursts (Motch et all 120071) . Thus, the 
most likely interpretation of the available information is that 
IGR JO 1363+66 10 is a Be X-ray binary with a neutron star 
(or perhaps a black hole as discussed below) accretor. In the 
following, we discuss IGR JO 1363+66 10 in this context. 



4.2. Accretion Regimes and the X-ray Luminosity 

Models for X-ray emission from Be X-ray binaries 
are based on the picture of a neutron star with a rel- 
atively strong magnetic field {B ~ 10 12 G) and a rela- 
tively slow rotation speed (P sp i n ~ 0.1-1,000 s) periodically 
accreting from the circumstellar disk around th e Be star 
(Rappapor t & van den Heuvelll982t ICorbetlll986l) . In most 
systems, the orbits have non-zero eccentricity, and "Type 1" 
X-ray outbursts occur when the neutron star makes its closest 
approach to the star so that the outburst periodicity is equal to 
the orbital period, P or b. 

Within this physical picture, there are two accretion regimes 
that are usually considered. One regime is at higher mass ac- 
cretion rates when the accretion pressure overcomes the cen- 
trifugal magnetic field barrier, and matter is accre ted directly 
onto the poles of the neutron star (Corbet 1996). The sec- 
ond regime is during times when the accretion is inhibited by 
the magnetic field. This is o ften called the propeller regime 
dlllarionov & Sunvaevl 19751) . and the X-ray emission is much 
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Table 3 

Spectral Results for XMMU J013549. 5+661243 



Parameter 




Value 


Power-law model 


N H 




(1.5^) X 10 22 cm" 2 


r 




1.4 ±0.3 


Fpt (unabsorbed, 0.2- 


-12keV) (1.5 


10- 13 ergs cm- 2 s-' 


X 2 /" 




23.1/26 


Blackbody model 


N ti 
kT 




(0.3;g |) x 10 22 cm- 2 
I.33ig-J| keV 






(3.3+J-f) x 10- 3 




21.5/26 


Bremsstrahlun; 


> model 


Nh 




(1.5!^) X 10 22 cm" 2 


kT 




>12keV 


^bremss " 




(2.9^ 5) x 10- 5 cm- 5 


x> 




' 22.6/26 



"The errors on the parameters are for Ax 2 = 2.7, corresponding to 90% confidence for one parameter of interest. 

*The normalization for this model is parameterized in terms of the radius of a spherical emitting region in units of km (R^ m ) and the distance in units of 10 kpc 

r The normalization for this model is 2.4 X I0~ l6 d 2 J n e n,dV cnr 5 where d is the source distance, and n e and n, are, respectively, the electron and ion number 
densities within the volume V. 



lower in this regime because the matter being accreted only 
reaches the neutron star magnetosphere rather than falling 
onto the neutron star surface. A gap is predicted between the 
minimum luminosity in which the system can be in the di- 
rect accretion regime and the maximum l uminosity in which 
the system can be in the propeller regime dCorbelll 19961) . and 
we define A to be equal to this luminosity ratio. For typical 
neutron star properties (1.4 Mq mass, 10 km radius, 10 12 G 
magnetic field), A = 170(P spin /l s) 2/3 dCampana et al.ll2002l) . 
and values of A r ange from ^1 00 for fast rotators to ~ 10,000 
for slow rotators (ICorbejl 19961) . 

Assuming a distance of 2 kpc and a spectral shape consis- 
tent with the power-law model measured using XMM-Newton, 
IGR JO 1363+66 10 has been observed by INTEGRAL at L x = 
1.04 x 10 35 ergs s" 1 in 2004, by Chandra at <1.4 x 10 31 
ergs s" 1 in 2007, and by XMM-Newton at 9.1 x 10 31 ergs s" 1 
in 2009. These luminosities are unabsorbed values measured 
in or extrapolated into the 0.2-12 keV band. The 2004 out- 
burst luminosity is at the lower end of typical values for nor- 
mal ("type I") Be X-ray bina ry outbursts, which are in the 
10 3S to 10" ergs s" 1 range (IStella. White & Rosned fl986t 
I Wilson et all 120051) . Still, it is likely (but perhaps not cer- 
tain if the neutron star spin or magnetic field strength have 
extreme values) that the source reached a high enough ac- 
cretion rate to ent er the direct accretion regime. Following 
ICorbetld 19961) and|Campana et al. (2002), we derive an upper 
limit on the luminosity produced by magnetospheric emission 
of L m < 1.04 x 10 35 /A = 6.1 x 10 32 (P spin /l s)" 2 / 3 . Thus, un- 
der the assumptions that the neutron star properties (mass, ra- 
dius, and magnetic field strength) are typical and the system 
reached the direct accretion regime during the 2004 outburst, 
the luminosity measured by XMM-Newton in 2009 implies a 
neutron star spin period less than 17 s for IGR JO 1363+66 10. 

While the propeller regime may extend to very low lumi- 
nosities, a third regime to consider is when the mass accre- 
tion rate onto the neutron star drops to zero. In this case, 
any X-ray emission produced could have contributions from 
the Be star and from the neutron star. Emission from the 



neutron star would be thermal in origin and would be ex- 
tremely soft. Such compo nents are seen at temperatures near 
50-150 eV in LMXBs ferown. Bildsten & Rufledgel [19981 
Camp ana et al.lfl998t iDegenaar et al.ll2009l) and although it is 
unclear whether this regime occurs for Be X-ray binaries, it 
has probably been seen for the HMXB Supergiant Fas t X-ray 
Transient (SFXT) IGR J17544-2619 din 'i Zand! 12005). Dur- 
ing the XMM-Newton observation of IGR JO 1363+66 10, we 
are not seeing this regime since we measure a hard spectrum, 
which is either non-thermal or has a much higher temperature 
(kT =1.3 keV) than would be expected from a non-accreting 
neutron star. Since we did not detect the source during the 
Chandra observation, the spectrum is not constrained. As- 
suming the same A^h as seen during outburst (A^h = 1 .6 x 10 22 
cm -2 ), we can put limits on the flux and the temperature of 
the neutron star by requiring an ACIS count rate below the 
upper limit given above. We use the magnetic neutron star 
model nsa (Pa vlov et al.lll995l) . assuming emission from the 
entire neutron star, a mass of 1.4M Q , a radius of 10 km, a 
distance of 2 kpc, and a magnetic field of B = 10 12 G. The 
upper limit on the (unredshifted) surface temperature is then 
80 eV, and the unabsorbed 0.2-12 keV luminosity upper limit 
is <2.4 x 10 32 ergs s" 1 . Changes in the assumptions will affect 
this, particularly changes in the assumed A^h or distance. For 
instance, a choice of A^h = 2 x 10 22 cm -2 gives kT < 85 eV 
and L x < 3.2 x 10 32 ergs s" 1 , or a distance of 3 kpc gives 
kT < 88 eV and L x < 3.7 x 10 32 ergs s -1 . If the X-ray emis- 
sion in this quiescent state is dominated by hot spots at the 
polar caps, their temperature may be higher, but the tempera- 
ture of the rest of the neutron star surface must be even lower. 

4.3. The Quiescent Luminosity and Possible Implications 

The above analysis suggests that the IGR JO 1363+66 10 
luminosities can be explained within the standard picture 
for Be-neutron star X-ray binaries, but the luminosity up- 
per limit that we infer from the Chandra observation is 
lower th an has been p re vious ly reported from these sys- 
tems. Camp ana et al.l d2002l) present a study that fo- 
cuses on quiescent X-ray observations of 3 Be X-ray bina- 
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ries, and their luminosities are (l-3)xl0 35 ergs s" 1 (0.1- 
10 keV) for A 0538-66, (0.8-2) x 10 33 ergs s" 1 (0.5-10 keV) 
for 4U 0115+63, and 5 x 10 32 ergs s~ f (0.5-10 keV) for 
V 0332+53. In addition, quiescent luminosities of (3- 
9)xl0 33 ergs s" 1 and (2-4.5) xlO 33 ergs s -1 have been re- 
ported for GRO J2058+42 and A 0535+ 26, respectively 
dWilson et all 120051: INegueruela et~ai1l2000l) . The candidate 
Be X-ray binary XTE J 1829-098 has be en observed at a lu- 
mino sity of 3 x 10 32 (c// 10 kpc) 2 ergs s" 1 (lHalpern & Gotthelfl 
120071) . which is low for a Be X-ray binary, but still not as 
low as we find for IGR JO 13 63+66 10. It should, ho wever, 
be noted that, for A 0535+26. lNegueruela et"aT1 (l2000h found 
the source at this luminosity after it had previously been re- 
ported t o have a quiescent luminosity two orders of magnitude 
higher dMotch et al.lll99lT) . Thus, it is clear that even when 
they are not in outburst, these sources exhibit a large amount 
of variability, and many observations may be required for each 
source to define a quiescent luminosity. 

While the Chandra upper limit of <1.4 x 10 31 (d/2 kpc) 2 
ergs s" 1 obtained with a hard spectrum as seen in outburst is 
significantly lower than has been found for other Be X-ray 
binaries, this conclusion does depend on the vali dity of the 
distanc e determination. The 2 kpc distance from lReig et al.l 
(2005) depends on the measurement of t he V-band magni- 
tude, which is given in Reig et all (120051) to high precision 
(V = 13.29 ±0.02, based o n three mea s ureme nts), the extinc- 
tion, and the spectral type. iReig et all (|2005) give a value of 
E(B-V) = 1.6, and our determination is consistent with this 
suggesting that the uncertainty in the distance due to the mea- 
surement of extinction is not large. Pr obably the largest uncer- 
tainty is related to the spectral type. Reig et al. (2005) use a 
spectral type of B1V and an absolute magnitude of My = -3.2 
for the 2 kpc estimate, but they indicate that the star could be 
a more luminous sub-giant (i.e., B1IV). The difference in the 
absolute magnitude of a main sequence B-typ e star comp ared 
to a B-type sub-giant is 0.6-0.7 magnitudes (Cox 20QO), so 
it is unlikely that the star is more luminous than My ~ -4.0, 
which corresponds to a distance upper limit near 3 kpc. This, 
in turn, indicates an upper limit on the quiescent X-ray lu- 
minosity of <3.2 x 10 ergs s -1 , which is still significantly 
lower than the X-ray luminosities given above for compari- 
son. 

Although it is not clear why the quiescent luminosity is 
so much lower for IGR JO 1363+66 10, there are several po- 
tentially interesting possibilities. Among the more trivial is 
the possibility that the short (5 ks) Chandra observation oc- 
curred during an X-ray eclipse. Although eclipses are not 
common for Be X-ray binaries due to their typically wide or- 
bits, we currently cannot rule out this possibility since the 
orbit and binary inclination for IGR JO 1363+66 10 are un- 
constrained. It is also possible that Chandra did not detect 
the source either because the compact object was highly ob- 
scured during the Chandra observation or that Chandra hap- 
pened to catch the source during a low point in its normal 
short-term variability. Finally, there is also the possibility 
that the circumstellar disk, which was strongly present in 
2004 and 2009 based on the detection of the Ha emission 
line with an EW of ^-50A, had dissipated during the 2007 
Chandra observation. Although dissipation and re-formation 
of a circumstellar disk on this time sc ale did occur in the 
case of GRO J 1008-57 dCoe et al.ll2007l) . and significant disk 
loss was o bserved for IGR J06074+2205 on a time scale of 
~3 years (IReig. Zezas & Gkouvelisll2010l) . in both of these 



cases, the Ha lines were significantly weaker than measured 
for IGR JO 1363+66 10. There is no precedent that we are 
aware of for dissipation and re-formation of a circumstellar 
disk in a system with an Ha EW as large as the one mea- 
sured for IGR JO 1363+66 10. While there are possible expla- 
nations for the luminosity upper limit measured over a short 
period of time by Chandra, these explanations do not apply 
to the XMM-Newton observation. The luminosity measured 
by XMM-Newton would also make it one of, if not the, lowest 
luminosity Be X-ray binaries. Thus, it is likely that the low 
luminosities measured for IGR JO 1363+66 10 require that this 
system has some unusual physical properties. 

The Be X-ray binaries with known compact object type 
have neutron stars, but binary evolution models predict 
that between zero and two of the 64 known Be X-ray bi- 
nary systems harbor a black ho le instead of a neutron star 
dBelczvnski & Ziolkowskil [20091) . The observational signa- 
tures for Be X-ray binaries with black holes are not as clear as 
neutron star signatures, such as pulsations. At high mass ac- 
cretion rates, one would expect the X-ray emission from black 
holes to be softer and brighter due to the presence of an inner 
accretion disk that cannot form in systems with highly mag- 
netized neutron stars due to the magnetosphere. At moderate 
to low mass accreti on rates, accreting black holes typically 
have hard spectra dRemillard & McClintockll2006l) . At very 
low mass accretion rates, black hole systems would likely 
be fainter than neutron star systems because the latter exhibit 
magnetospheric and surface emission, which are both absent 
from black holes. The expected hard spectrum and the faint 
quiescent luminosity are both consistent with the observed 
properties of IGR JO 1363+66 10. 

There may, however, be a less exotic explanation if the 
binary orbit of IGR JO 1363+66 10 causes a neutron star to 
sample an unusually low density part of the circumstellar 
disk. This could occur if the system has a long or bital period 
with re latively low eccentricity. As pointed out by lReig et al.l 
(2005), the EW of the Ha line for IGR JO 1363+66 10 is 
one of the strongest of any Be X-ray binary, suggesting the 
presence of a large circumstellar disk, which might favor a 
large and nearly circular orbit over a highly eccentric or- 
bit. While the system may h ave a large circumstellar disk, 
INegueruela & Okazakil ([2001) have shown that tidal itera- 
tions of the neutron star can lead to truncation of the circum- 
stellar disk in low eccentricity systems. In addition to provid- 
ing a possible explanation for the low quiescent luminosity of 
IGR JO 1363+66 10, the fact that the disk can be truncated to a 
size smaller than the Be star's Roche lobe means that the sys- 
tem w ill not show luminous outbursts dOkazaki & Negueruelal 
1200 lb . which is consistent with what has been seen so far for 
IGR JO 1363+66 10. 

Like IGR JO 1363+66 10, the orbital parameters for the 
other Be X-ray bi naries with the strongest Ha emission 
lines, A 1118-61 6 dCoe et all 1 1994b and IGR J01583+6713 
dKaur et al.ll2008l) . are not known. It has been suggested that 
A 1 118-616 has a nearly circular orbit based on it s having a 
small number of outbursts with no clear periodicity dCoe et alj 
1994). It has been argued that both of these systems have rel- 
atively long orbital periods (^0.6-2.2 years) based on mea- 
sured neutron star spin periods of 405.6 s and 469.2 s and 
the f Q rb-fi P in rel a tionship for Be X - ray binaries dCorbetl 19861 : 
lCoeetalJll99lt iKaureta l. 2008; iDoroshenko et all 120101) . 
However, for IGR J01583+6713, the detection of X-ray pul- 
sations was noted as being marginal in lKaur e t al. (200g). Al- 
though there is no direct indication of the orbital or spin pe- 
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riod for IGR JO 1363+66 10, it should be noted that the limit of 
^spin < 17 s mentioned above could be taken as evidence for a 
smaller P or b; however, we view the evidence for a limit on the 
spin period as being relatively weak. 

5. SUMMARY AND CONCLUSIONS 

In summary, the detection of the variable hard X-ray source 
XMMU J013549.5+661243 at the location of a Be star pro- 
vides confirmation that IGR JO 1363+66 10 is a Be X-ray bi- 
nary. Although some of its properties suggest that it is a fairly 
typical HMXB, its 2007 non-detection by Chandra indicate a 
quiescent X-ray luminosity that is significantly lower than has 
been measured previously for other Be X-ray binaries. While 
some possible explanations for the low luminosity (an eclipse, 
a large change in the column density or spectrum, an extended 
drop in mass accretion rate, or dissipation of the Be star's 
circumstellar disk) may be consistent with relatively normal 
Be X-ray binary properties, other possibilities would require 
that IGR J01363+6610 has unusual properties. One possibil- 
ity is that the system has a large and relatively circular orbit, 
which could explain the low duty cycle for outbursts and, per- 
haps, the low quiescent luminosity. Another very interesting 



possibility is that this could be a Be-black hole system. The 
former suggestion could be confirmed by a measurement of 
the orbital period (although this will be challenging since the 
source does not seem to produce regular outbursts), and the 
latter suggestion could be refuted with the detection of pulsa- 
tions during another outburst from the source. 
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